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USE OF REMOTE SENSING FOR LAND USE POLICY FORMULATION 


Introduction 

This has been an exciting year for the Center with many significant accom- 
plishments. Also, with the evolution of the new global habitability programs, 
particularly the land-related sciences efforts, we have found that many of our 
existing programs fit well into this new proposed thrust area of NASA. In addi- 
tion, preliminary analysis of thematic mapper data shows that it will open broad 
new horizons in providing informational content that could not before be 
obtained from available remotely sensed data. 

Because of the diversity of the Center's activities and the wide nunber of 
participating scientists, our research and technology transfer efforts involve a 
number of subject matter areas. As renewable resources become increasingly 
important, the need for identifying forest types becomes extremely important. 
The study by Carl Ramm and Denny Hudson examined this area with the use of 
Landsat data. Their overall accuracy for forest cover type identification was 
85H. Thus, these and future techniques will be extremely useful to many indus- 
tries. Further research efforts will study forest vegetation identification as 
it relates to biomass, deforestation and NCOg flux studies. 

Energy balance and hydrologic cycles require knowledge of synoptic tempera- 
ture patterns. Extensive analysis showed that HCMM data can provide this type 
of information and it appears that GOES data would also provide valuable infor- 
mation on a diurnal basis. Preliminary indications are that this information 
will help with the understanding of man's impact on the surface conditions of 
this planet and the data can be used in models to aid in determining energy bal- 
ance and hydrologic studies. 

As man's activities (with increasing population demands throughout the 
world) are impacting land cover and land use, it is increasingly important that 
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we can develop techniques for detecting change and updating inventories of sur- 
face conditions. The project entitled Land Surface Change Detection outlines 
significant progress in this general and important area. This analysis has 
already shown classification accuracies for the urban and forestry test sites 
approaching the 75 % accuracy standard used in the preliminary study. A much 
better understanding is evolving relative to the relationship with the diversity 
and spatial complexity of the test sites being classified. The preliminary 
thematic mapper data available shows that significant increases in accuracy will 
be possible with this information. 

Soil mapping is one of the major undertakings in inventorying our key 
resources worldwide. The study on "Determination of Soil Hap Unit Composition 
by Electronic Scanning Densitometry" examined new methods for aiding in soil map 
unit composition determinations and outlines the accuracy that can be expected. 

The research dealing with virus infection of vineyards and blueberry fields 
is an important basic study. Plant diseases have a major impact on crop produc- 
tivity and fundamental studies are required, such as this one, to better del- 
ineate the value of remote sensing, its limitations and research directions for 
improving sensor capabilities. 

Several contractual activities were done under the broad responsibilities 
of the Center for Remote Sensing. These included an Important Farmlands Inven- 
tory for several Michigan counties which identified areas of prime and unique 
farmlands. Another project studied changes in the aquatic vegetation in Saginaw 
Bay and showed how remote sensing could be a key tool in delineating over time 
the extreme variations that occurred in aquatic vegetation areas. The study 
used historic aerial photography and the capabilities of a geographic informa- 
tion system to analyze a 28-year time sequence of changes in the amount and dis- 
tribution of aquatic vegetation. The analysis showed clearly that dramatic 
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shifts do occur and it provided aerial statistics and documented some of the 
cyclic patterns involved. A special study was undertaken using aerial photogra- 
phy in hybrid poplar research. The hybrid poplar research program is particu- 
larly significant as we look to the future when bioenergy plantations will 
become increasingly important. These efforts are helping to develop techniques 
for inventorying the bio-resources available through hybrid poplar plantations. 

Close ties have evolved with an international project in the Department of 
Resuorce Development at MSU which is titled Comprehensive Resource Inventory and 
Evaluation System (CRIES). Major resource inventory projects are on-going in 
Jamaica, Honduras, Dominican Republic and Kenya. This project provides support 
to the Center for certain CRIES-related activities, which have many implications 
for the global habitability program from an international standpoint. Through 
the CRIES project new detailed data bases have been developed for several coun- 
tries and cross-checked with ground truth information. Our past international 
work could provide key sites for many studies on changes in cropping practices, 
deforestation, and in the broad areas of biological productivity and impacts on 


the hydrological cycle 
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RESEARCH PROJECTS 



IDENTIFICATION OF CONIFEROUS FOREST TYPES IN MICHIGAN 
USING LANDS AT IMAGERY 

Carl W. Ramm 
Department of Forestry 

William D. Hudson 
Center for Remote Sensing 
Department of Forestry 

This study evaluated the use of Landsat computer enhanced Imagery for map- 
ping coniferous forest types in Michigan's northern Lower Peninsula. Visual 
interpretation procedures were developed and tested over two sites to determine 
the feasibility of Identifying coniferous species. The accuracy achieved by two 
interpreters was compared and sunmarized in contingency tables. Overall clas- 
sification accuracies were 85 and 73 percent whereas individual species 
interpretability accuracies ranged from a low of 32 percent for mixed pine 
stands ;o 95 percent for jack pine plantations. Most of the errors in mapping 
the pines were accounted for by confusion between the individual pine species. 
The swamp conifer type had consistently low interpretation accuracies at both 
test sites. Other factors affecting interpretation accuracies were also identi- 
fied and are reported. 

The objective of this study was to evaluate the accuracy of mapping coni- 
ferous forest types in the northern Lower Peninsula of Michigan through visual 
Interpretation of computer-enhanced Landsat imagery. 

The climate, physiography and soils of the northern Lower Peninsula combine 
to create a unique habitat region, as compared to the southern Lower Peninsula 
or the Upper Peninsula. One of the most distinctive landscapes are high outwash 
plains near the center part of the region, which are dominated by stands of 
pines and oaks. These extensive, flat, sandy plains, which formerly supported 
the large white pine forests of Michigan, support much of the nearly seven mll- 
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lion acres of forest found In this unit* 

Hardwood forests are more extensive than the softwoods with aspen-birch, a 
sub-climax type, the most common throughout the region. The northern hardwood 
type Is most prevalent in the northwestern counties, while lowland hardwooda are 
primarily restricted to riparian sites. The softwoods, or oonlferous forest 
types, occupy approximately 22 percent of the region. Three-fifths of these 
conifers are pines; Jack pine, the most important pine sub-type, acoounts for 
more than half the acreage. Jack pine occurs as a relatively pure type in a 
broad belt from the central northeastern area, south and west to the southwest 
central area. Significant areas are also found in the east and westcentral 
counties. The swamp conifer types, consisting predominately of cedar with 
lesser areas of black spruce, balsam fir-white spruce end tamarack, comprise 
under 9 percent of the forest land. These swamp conifers occur most frequently 
as small patches of a few acres on wet lands. The exceptions are in the 
northeastern counties where these species are more frequent and exist as larger 
stands. 

Two test sites were chosen in the northern Lower Peninsula to be represen- 
tative of areas now supporting large acreages of conifers. The first test site 
was located in west central Wexford County (T. 22 and 23 N., R. 12 W.) and is 
underlain primarily by stratified sand and gravel outwash deposits confined to a 
broad valley (i.e., a valley train). Host of the area was cleared for agricul- 
ture at one time, but later was abandoned as unsuitable for sustained crop pro- 
duction. The present forest covers over 40 percent of the area and is predom- 
inately pine plantations. Red pine accounts for 60 percent of the plantations, 
jack pine is predominate in 29 percent of the plantations and 11 percent are 
mixtures of red and jack pine. Ten percent of the forest land is composed of 
the swamp conifer type. These stands, composed of scattered northern white- 
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cedar intermixed with lowland hardwoods , are concentrated along several creeks 
traversing the region. The second teat site waa located in northeastern Craw- 
ford County and Southeastern Otsego County (T. 28 and 29 N., R. 1 W.). This 
area is part of an extensive outwaah plain and is typical of the "jack pine 
flats" of the northern Lower Peninsula of Michigan. No evidence of land clear- 
ing for agriculture is present and the forest is entirely natural (e.g. there 
are no plantations). The area is nearly entirely forested* with conifers dom- 
inant on 72 percent of the area, the remainder being predominately hardwoods, 
grass or brush. Pines are located throughout the site except for two large 
swamps which support lowland conifer species. Jack pine is by far the predom- 
inant species throughout the site and represents over 70 percent of the softwood 
acreage. 

Since the study dealt only with coniferous vegetation, a "leaf-off" scene 
was chosen because of its sharp tonal (color) contrasts of several coniferous 
forest types. Based upon an initial analysis of this scene and the characteris- 
tics of the coniferous forest distribution in the northern Lower Peninsula, a 
classification scheme was developed. 

Because of their large aggregate acreage and wide distribution, stands of 
Jack pine, red pine and mixed stands were delineated but white pine, which 
represents less than 1 percent of the forest land, was not separated out. Low- 
land conifers, because of their small total acreage and tendency to form highly 
mixed stands, were grouped into a single category. 

To facilitate interpretation and mapping functions, a pin-registered over- 
lay system was utilized. For each test site a base map was constructed from 
U.S. Geological Survey topographic maps. These maps, at a scale of approxi- 
mately one inch to the mile, were registered, via punched holes aligned by a 
metal pin bar, with the "ground truth" maps and the Landsat interpretation maps. 



This system permitted prsclst ovsrlaying and comparison of two or nor* of tha 
separata maps. 

Prior to actual photo-interpretation, tha two interpreters were given prel- 
iminary training, including the development of photo-interpretation training 
aids. Photo keys were prepared to illustrate the appearance of the different 
coniferous forest types on Landsat false-eolor composites. Additional training 
consisted of the systematic comparison of several examples of each forest type 
on high-altitude color Infrared photography and on the Undaat color composite, 
coupled with specifically gathered forest stand Inventory measurements. 

The Interpretation entailed the enlargement of the Landsat color composite 
on a precision rear projector. Magnification was controlled by matching the 
enlarged Landsat scene to the previously-prepared base map, and then replacing 
the base map with a blank sheet of polyester film. Utilizing the previous photo 
comparisons and photo keys, the Interpreters identified and delineated the boun- 
daries of all coniferous forest stands. Tone, end, in some cases, texture, were 
tne mein Interpretation criteria for the pine classes. Proximity to stream 
courses and water bodies were an additional aid in distinguishing between upland 
pine types and lowland swamp conifers. 

To evaluate the accuracy of the visual interpretation procedure, the 
Landsat Interpretations were compared with previously compiled cover type maps. 
These maps, tdilch were prepared specifically for this project, were constructed 
from photo-interpretation of medlim-scale (1:24,000) color infrared photography. 
Additionally, both U.S.D.A. Forest Service and Michigan Department of Natural 
Resources forest cover type maps were consulted in conjunction with ground 
verification by field crews. 

All errors in each Landsat interpretation map were identified by superim- 
posing these maps on the cover type maps using the pln-reglstratlon system. The 
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*rea of eaoh stand was maasurad using a 160 dot/lnot)* grid and the raaults sun- 
narltad In eontlnganoy tablas (Tablas 1 and 2). 


Tabla U Landsat Clasalfloatlon Performance, 
~ Waxford County Tast Slta. 


Known 

Covar 

Typa 

Number 

Rad 

Pina 

i 

of Sample Points Classified as— 

Jaok Pina Swamp Non- 

Pine Mixtures Conifers Forest 

Covar 

Typa 

Total 

Rad pine 

m 

8 

130 

26 

103 

1,041 

Jack pine 

31 

121 

17 

10 

«8 

497 

Pir.o mixtures 

71 

7 

21 

9 

13 

187 

Swamp conifers 

16 

2 

0 

m 

36 

197 

Non-forest 

57 

3 

36 

89 

2,590 

2,775 

Total 

"553 

Ml 

276 

S77 

2,790 

"TST 


Percent correct* 

82 

95 

32 

52 

93 

85 3 


Values along tha diagonal raprasant eorraetly dallnaatad and ldantlflad oovar 
typaa. 

2 

Ratio of diagonal valua to tha total count of that oovar typa as infcerpreteo 
from tha Landsat seana. 

^Overall classification accuracy; ratio of tha sun of diagonal valuas to tha 
total nunbar of saupla points. 

Thasa two-way crosr»tabulatlons paralt tha rasulta to ba vlawad fron two 
approachas. By reading across a row ona may obsarva vhlch catagorlas, and to 
what dagraa of accuracy t a particular covar typa was dallnaatad. By raadlng 
down a column, ona may obsarva tha actual covar typaa composing aach Landsat 
classification. For discussion purposes, accuracy was summarized as the ratio 
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(expressed as a percent) of correctly delineated and classified stands to the 
total count of that cover type as interpreted from the Landsat soene. 


Table 2. Landsat Classification Performance, 
Crawford County Test Site. 


Number of Sample Points Classified as— 1 
Known Cover 

Cover Red Jack Pine Swamp Non- Type 

Type Pine Pine Mixtures Conifers Forest Total 


Red pine 

23 

38 

0 

* 

« 

16 

78 

Jack pine 

9 

1,500 

18 

18 

33 

1,578 

Pine mixtures 

0 

23 

33 

11 

10 

77 

Swamp conifers 

0 

125 

1 

398 

19 

543 

Non-forest 

2 

301 

1 

22^ 

07 

833 


Total 

34 

1,987 

53 

650 

385 

3,109 


2 

Fercent correct 

68 

75 

62 

61 

80 

73 3 


^Values along the diagonal represent correctly delineated and identified cover 
types. 

2 

Ratio of diagonal value to the total count of that cover type as interpreted 
from the Landsat scene. 

^Overall classification accuracy; ratio of the sum of diagonal values to the 
total number of sample points. 

Table 1 sunmarlzes the Landsat interpretation accuracies for the Wexford 
County test site, which had an overall classification accuracy of 85 percent. 
The high accuracies obtained for both red and jack pine are probably a direct 
result of their hifch stocking levels and frequent occurrence in pure planta- 
tions. Mixed plantations, where neither species represents more than 75 percent 
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of the stocking, were consistently mlsclasslfled. Further Investigation 
revealed that two-thirds of these errors were accounted for by commission and 
omission errors to the red pine type. Additional cross checks confirmed that 
the majority of errors in the pine types were due to misclasslflcations among 
the several pine species. 

Accuracies obtained without considering the wlthln-pine variability, that 
is, utilizing a "pooled" pine category, show that pines were correctly inter- 
preted 93 percent of the time. 

The swamp conifer type was correctly delineated and classified with approx- 
imately 73 percent accuracy. The lower reported accuracy (52 percent) resulted 
from a large commission type error. Upon examining the medium-scale photogra- 
phy, the majority of these errors were found to occur over swamp areas dominated 
by lowland brush and/or lowland hardwood species. Further "lumping" of species 
was attempted and indicated that forest land as a single category was inter- 
preted with 90 percent accuracy. 

Table 2 summarizes the Landsat interpretation accuracies for the Crawford 
County test site. The lower accuracies for both red and Jack pine may be caused 
by their lower stocking levels, with many stands grading into non-forest 
categories (i.e., less than 25 percent stocked). The various forest typos 
encountered at this site exhibited less contrast than those at the Wexford 
County site. Most errors in the pine categories occured among the various pine 
types, and in recognizing hardwoods intermixed with low density pines such as 
jack pine. Again, "pooling" the pine types increased overall interpretation 
accuracy to 77 percent for the pine category. 

Errors in swamp conifer type were similar to those found in the Wexford 
County site. Lowland brush, often with scattered trees, were the most fre- 
quently misinterpreted areas, although some confusion between jack pine and 
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swamp conifers did occur. The overall forest/non-forest accuracy was 81 per- 
cent. This was substantially below that obtained over the Wexford County site. 




INVESTIGATION OF SYNOPTIC TEMPERATURE PATTERNS IN MICHIGAN AS 
DETERMINED VIA VIES AND HCMM THERMAL IMAGERY 

Jon Barthollc 
Center for Remote Sensing 
Department of Resource Development 

Stuart Gage 

Department of Entomology 

Ardeshir Goshtasby 
Department of Computer Science 

Considerable additional progress had been made in the use of meteorological 
satellite data information. Additional digital GOES data tapes have been 
acquired. These took a considerable length of time to obtain, however, they 
have now been examined, screened and are ready for detailed statistical 
analysis. 

HCMM data has been acquired of approximately 12 scenes over Michigan, one 
scene over Georgia and two scenes over Florida. Preliminary analysis of this 
information shows that the data looks good and provides extremely valuable 
information on temperature patterns and reflectances as influenced by man's 
activities and alterations of the earth's surface. 

Use of this information in more quantitative approaches has resulted in a 
paper (see Appendix). This paper pulls together past unpublished information by 
the author as well as new analysis that has been made during the last contract 
period. The paper sunmarizes the past information, our new knowledge and the 
implications and neei' for further analysis. A part of this activity led to a 
proposal entitled, "Temperature and Reflectance Monitored from Satellites as an 
Indication of Shift and Impact of Vegetation Change." This proposal will be 
funded out of NASA Goddard. 


13 


LAND SURFACE CHANGE DETECTION USING 


i 


SATELLITE DATA AND A GEOGRAPHIC DATA BASE 

William R. Enalin 
Center for Remote Sensing 

Anil Jain, Eric Backer and Ardeahir Goahtasby 
Department of Computer Science 

Richard Hill-Rowley 
Department of Geography 

Michael Scieszka 

Michigan Department of Natural Resources 

The basic objective of this research is to develop methods and techniques 
for the detection, monitoring and modeling of land surface changes from data 
acquired by satellite sensor systems. Data from the Landsat multispectral 
scanner (MSS), weather satellites, and the new Landsat Thematic Mapper (TM) are 
being evaluated to determine the capability to accurately identify land cover 
and land use types and detect land surface changes. 

The research is being conducted in cooperation with the Michigan Resource 
Inventory Program, which is currently implementing a legislated statewide infor- 
mation system. The system Includes digitized photo-derived land use information 
which will serve as the baseline data set for comparisons with computer-aided 
classifications of satellite digital images. The aim is to develop and evaluate 
techniques and methodologies for integrating and processing data obtained from 
satellite sensors with data layers in the geographic information system of the 
state. This research will eventually lead to the design of multi-stage change 
detection research. 

Several research areas have been identified which will provide techniques 
needed as components of an overall system for monitoring land cover/use changes 
and updating regional and statewide geo-referenced data bases. Other areas of 
research are being defined related to integrating and analyzing satellite data 
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with selected earth surface attribute sets of the data base, particularly infor- 
mation on forest stand type, size and density. A summary of the current efforts 
in each research area follows. 

IMAGE REGISTRATION 

Accurate image registration is a prerequisite for most change detection and 
inventory update methods. We have, therefore, concentrated our initial work in 
this area. Techniques to handle three types of image ‘ registration will be 
required for the research, i.e., image to image, image to geographic data base 
and image with ancillary data sets such as digital elevation data. 

Several multi-image registration algorithms have been developed and 
evaluated during this reporting period. Emphasis had been placed on automatic 
digital image registration procedures that handle images with translational, 
rotational and scaling differences. One procedure, which has been reported in a 
paper (see Appendix), involves 1) segmentation of images, 2) isolation of dom- 
inant objects, 3) determination of corresponding objects in the two Images and 
4) estimation of transformation parameters using the center of gravities of 
objects as control points. This method has been used to register different 
Images of Michigan acquired by the Heat Capacity Mapping Mission (HCMM) satel- 
lite. HCMM images have also been registered with Images from the Geostationary 
Operational Environmental Satellite (GOES) using invariant moments. 

Attempts were made to register airborne thematic mapper (TM) simulation 
data with a Landsat MSS scene. Geometric distortions in the aircraft data 
prevented a good fit, however we believe that a recently-acquired Landsat TM 
scene can be properly registered to Landsat MSS scenes. 

We have also begun to explore computer-based techniques for the determina- 
tion of pixel elevation based upon the registration of stereo Images. These 
techniques will be useful in the analysis and classification of data from 
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polntable sensors on satellites. 

Work on the registration of digital images to the states geographio infor- 
mation system (GIS) data files has not been started due to problems in reformat- 
ing the GIS data for tape output. When tape reformatlng is accomplished y rou- 
tines will be written to handle polygon-to-raster conversions and image to data 
base registration. 

As an initial check of registration potential , a county area on a Land sat 
TM Band 5 scene was photographically enlarged to a scale of 1:48,000 and over- 
layed with a clear-acetate copy of a computer-drawn land use map from the Michi- 
gan Resource Inventory System. The information content (tonal areas) in the 
Landsat TM scene coincided perfectly with the polygonal areas of the land use 
categories shown on the map. Based on the goodness of fit of the photographic 
products, we believe that the corresponding digital data sets can be easily and 
accurately registered. 

To date no work has begun on merging ancillary data sets with digital 
images from satellite sensors. However, digital elevation data for two 
1 : 250, 000 topographic quadrangles has been ordered and received from the 
National Cartographic Information Center. 

IMAGE SEGMENTATION 

Change detection through scene analysis of multiple digital images often 
requires adequate segmentation of the images. Various techniques for image seg- 
mentation will eventually be evaluated including edge operators, thresholding, 
region growing and multi- image segmentation procedures for determining 
corresponding regions in two data sets. Resultant binary masks for forest and 
urban features will then be employed to limit image classification and change 
identification within predefined areas. Preliminary research was conducted this 
reporting period on using image segmentation techniques to Isolate road networks 
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in aircraft thematic mapper simulation data. Roadways are important recognition 
criteria in the identification of residential subdivisions which are an impor- 
tant land use category requiring frequent updating. 

IMAGE CLASSIFICATION 

The Initial work in this research area is reported in a Ph.D. Dissertation 
entitled An Evaluation of Digital Landsat Classification Procedures for Land Use 
Inventory in Michigan . 1 

The major research objective of this study was to evaluate the informa- 
tional value of Landsat digital data in the context of providing land use data 
to users in Michigan. Information on user needs was developed through an exten- 
sive survey focusing primarily on land use planners. Land categories were iden- 
tified and subsequently modified in order to be compatible with processing of 
Landsat data. 

Three test sites, each characterized by a distinct land use type (agricul- 
tural, urban and forest) were chosen to evaluate Landsat performance. For each 
of these sites Landsat data for August 16, 1979 were classified by means of 
three commonly available algorithms (maximum likelihood, minimum distance-to- 
means, grouping of cluster classes) . Accuracy evaluation of the resulting clas- 
sifications included examining the effects of generalization for geographic 
information system formats. 

Conclusions derived from the study can be summarized in three major areas. 
Analysis of land use categories identified by questionnaire responses, in con- 
junction with limitations imposed by the characteristics of the Landsat data, 
show that the minimum land use requirements of the majority of planners surveyed 

1 Hill-Rowley, Richard, 1982. An Evaluation of Digital Landsat Classification 
Procedures for Land Use Inventory in Michigan . Ph.D. Dissertation, bepariment 
of Geography, Michigan State University, East Lansing, Michigan. 
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could not be fulfilled by the classification procedures. The categories that 
remained, nevertheless, represent the dominant land use classes In each of the 
areas selected. Accurate updating and/or remapping of these categories would be 
of considerable value. 

Classification accuracies for the urban and forest test sites approached 
the 75* accuracy standard used in the study. Accuracies for the agricultural 
site were significantly lower. Variations were present among algorithms, with 
the un supervised grouping of cluster classes substantially superior in the agri- 
cultural site, maximum likelihood the most accurate in the urban site and 
minimum distance-to-means preferable in the forest site. Despite this distribu- 
tion, no significant difference existed at the 95£-confidence level between the 
overall performance of the algorithms. 

Accuracy analysis of the generalized Landsat classes, when compared with 
geocoded land use types derived from airphoto interpretation, resulted in a sub- 
stantially lower accuracy than the same comparison between actual aerial photog- 
raphy and the pixel-based Landsat classification. The extent of this accuracy 
loss seems to have a direct relationship with the diversity and spatial complex- 
ity of the test site being classified. 

CHANGE DETECTION 

No work on change detection techniques was performed with Landsat digital 
data this reporting period due to the previously-stated problems with the tape 
routines of the state's GIS system. However, preliminary evaluations are being 
made of the superimposed Landsat TM/GIS Landuse Map product. The landuse data 
was obtained from 1978 color infrared aerial photography whereas the TM data was 
acquired in October 1982. The most predominant changes identified through 
visual analysis were forest clear-cut areas and new pads for oll/gas wells. The 
oil/gas pads are easily identifiable on TM images because of their regular ' 
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geometry and aeral dimensions. Work is underway to automatically pin-point the 
geographic location (latitude and longitude) of oil/gas pads. This will require 
an evaluation of the relative geometric aoouraoy of the state's GIS data files 
and the Centers Landsat geometric correction routines. 




DETERMINATION OF SOIL MAP UNIT COMPOSITION 


BY ELECTRONIC SCANNING DENSITOMETRY 
Delbert L. Mo km a 

Department of Chop and Soil Science* 

Modern aoll survey reports provide information on the soil resources of an 
area to a wide spectrum of public and private users. Natural soil boundaries 
are located by soil scientists in the field and delineated on an aerial photo- 
graphic base. Some map units are made up of primarily one kind of soil, while 
others are made up of two or more soil series. Most soil map units include 
small areas of other soils because the map scale and nature of the soils in the 
area do not permit the separation of these small areas. In some map units these 
Included soils have properties that differ substantially form those of the dom- 
inant soil or soils and, thus, could significantly affect use and management of 
the map unit. A soil complex is used to represent areas where two or more soils 
make up a large portion of the map unit but are so intricately mixed or so small 
in size that they cannot be separated on the soil map. 

The percent composition of each soil in a soil complex is reported in the 
soil survey report. These data are obtained from field observations. Statisti- 
cal sampling methods and ground enuneratlon techniques have been developed to 
improve the accuracy of these map uilt composition data. However, most of these 
methods are labor-intensive and not extensively utilized. 

Many soil complexes in Michigan have a mottled appearance on aerial photog- 
raphy because of differences in soil moisture, organic matter content, surface 
texture and/or vegetative variation. Imagery of bare or sparsely vegetated 
ground collected during spring months generally provides the most distinctive 
patterns. Soil scientists recognize relationships between some tonal patterns 
and kinds of soil. Thus, photointerpretation is an important part of soil map- 
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ping. Electronic scanning dsnsltonstry can olasslfy tha range of density values 
into range categories and rapidly determine the percentage composition within an 

area. 

The objective of this study is to evaluate the suitability of electronic 
scanning densitometry to accurately determine soil complex composition In a 
timely and cost-effective manner. 

METHODOLOGY 

Three soil complexes. Tappan-Londo loams. 0 to 2 percent slopes. Tappan- 
Avoca complex. 0 to 3 percent slopes and Guelph-Londo loams. 0 to 6 percent 
slopes, were selected for investigation. The classification of the soil series 

are: 

Tappan - Typlc Haplaquolls, flne-loany, mixed (calcareous), meslc 
Londo - Aerie Glossaqualfs, fine-loamy, mixed, mesio 
Avoca - Entic Haplaquods, sandy over loamy, mixed, meslc 
Guelph - Glossoborlc Hapludalfs, fine-loamy, mixed, meslc 
Two fields of each complex were located in Tuscola County, Michigan uhere a 
modern soil survey is In progress. Thus, the findings of this study would bene- 
fit the survey party in preparing composition estimates of these soil complexes. 
Each of the fields had no or sparse vegetation at the time the aerial photograph 
was taken. 

Using black-and-white panchromatic aerial photo transparencies the six 
fields were digitized using an EYECGM Image Scanner (Spatial Data Systems). The 
cell size used was about 3 m square. 

To determine the composition of the soil complexes, observations were made 
in a grid pattern at 30 m intervals. Soils were classified to the soil series 
level in the field. Surface horizon texture, type of surface horizon, depth to 
free carbonates, and thickness of sand, where present, were recorded for each 
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observation. 

The tolls data will ba overlaid on tha digitized density naps to datarnlna 
tha range of density values for each soil in the eonplex. Then the oonpoaltlon 
of the complexes will be determined using the digitized data and oonpared with 
that determined from field observations. 

RESULTS 

An example of the digitized soil cap for the Tappan-Londo eomplex is shown 
in Figure 1. The density ranged from 0 (dark) to 255 (light). The range of 
density for Tappan was 50 to 100 and for Londo was 100 to 174. For the field in 
Figure 1 the composition is 52 percent Tappan and 48 peroent Londo. 

The composition of the three soil complexes based on field observation is 
given in Table 1. The Tappan-Londo complex is composed of about 50 percent Tap- 
pan and 40 percent Londo with 5 percent Parkhlll and 4 percent other soils. 
Parkhlll had properties similar to Tappan while the other soils were more simi- 
lar to Londo. Combining these soils the composition of soils with reflectance 
properties similar to Tappan is 55 percent and those similar to Londo is 44 per- 
cent. These values are not significantly different from those estimated fcr the 
field in Figure 1 using electronic acenning densitometry. This indicates the 
potential for the use of this technique to determine composition of soil com- 
plexes. 

FUTURE WORK 

The comparison of the soils data with the density maps of the remaining 
five fields will be completed. This will give a more accurate evaluation of the 
use of electronic scanning densitometry to determine composition of soil com- 
plexes. 
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Table 1. Composition of three soil complexes as determined from field 
observations. 


Complex 


Soil Series* 


Percent 


Tapp an -Lon do Tappan 50 

Londo 40 
Parkhill 5 
Capac 2 
Other soils 2 

Tappan-Avoca Tappan 22 

Avoca 1 
Selfridge 48 
Wixom 2 
Capac 5 
Londo 3 
Parkhill 1 
Brookston 1 
Met amor a 1 
Corunna 3 
Belleville 8 
Other soils 5 

Guelph-Londo Guelph 52 

Londo 7 
Marlette 13 
Capac 9 
Me tea 6 
Owosso 3 
Metamora 1 
Sel fridge 3 
Oshtemo 2 
Parkhill 1 
Other soils 3 


*0ther soils represents soil series which were encountered at only one 
observation. 


DELIMITING AREAS OF VIRUS INFECTION IN VINEYARDS AND BLUEBERRY 
FIELDS OF SOUTHWESTERN AND WESTERN MICHIGAN 

Donald C. Ramsdell 

Department of Botany and Plant Pathology 

Adele M. Childress-Roberts 
Department of Botany and Plant Pathology 

Michigan blueberry growers have experienced considerable economic loss from 
the infection and spread of two serious virus diseases, blueberry shoestring 
virus (BBSSV) and blueberry leaf mottle virus (BBLMV). Approximately 56% of the 
commercial plantings are of the cultivar Jersey; one of the most susceptible 
varieties to both diseases. 

Another virus disease on grape, peach rosette mosaic virus (PRMV), has cost 
growers millions of dollars. More than 50 vineyards in the southwest section of 
the state experience crop loss due to this disease. All three diseases are 
transmitted readily by insect or nematode vectors. 

An increase in the number of fields becoming infected and resultant 
economic implications, makes early detection important. Control strategies for 
these diseases include rouging and replacement of bushes and vines exhibiting 
visual symptoms. However, a four year latent period is necessary before symp- 
toms are expressed, while infected plants serve as inoculum reservoirs. The 
ELISA (Enzyme- linked Immunosorbent Assay) technique while providing an accurate 
procedure for detecting individual bushes, is costly and impractical on a com- 
mercial basis. Therefore, the need exists for the development of a cost- 
effective, accurate and rapid method to discriminate between diseased and 
healthy plants. One approach is the use of remote sensing to observe spectral 
and geometric changes in plant canopies, using photography and field spectropho- 
tometric readings at various phonological stages. If practical, these methods 
would enable growers to rapidly assess the extent of infection before control 
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measures are initiated 


MATERIALS AND METHODS 
Ground Truth 

To correlate changes observed on photographic imagery or spectrophotometric 
readings, areas within infected blueberry fields and vineyards were selected 
based on a visual assessment of diseased and healthy conditions. Maps (Fig. 1) 
of these areas were constructed and the verification of infected plants was made 
using ELISA results. 

Aerial Photography 

Aerial photographs were taken of the previously mapped areas in southwest 
and western Michigan at various phenological stages during the 1981 and 1982 
growing season. Color, black/white IR and color infrared (CIR) transparencies 
(Kodak Aerochrome Infrared Film 2424 and 2443) were obtained using two ELM/500 
Hasselblad camera systems fastened to a belly-mount of a Cessna Skyhawk. 
Imagery was acquired at approximately 305 m (AGL) using a 70 mm format with 80 
mm lenses and Wratten 12 and 25 filters. 

Low Level Imagery 

Tower Photography 

To assess whether infected plants could be detected at lower altitudes, a 
35 mm camera was mounted from a 15 ft. hydraulic tower. Color-, CIR and mul- 
tispectral (70mm) film was taken of adjacent known diseases and healthy plants 
in the field. Vertical and angled photographs were taken at various exposure 
settings during full bloom, petal fall, fruit set and during senescence in 1981. 
Color and CIR negatives were analyzed with a color densitometer to detect den- 
sity differences between plants at various wavelengths. 
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Spectral Signature 

Changes in the raesophyll layer of the leaf, the primary site of infrared 
reflectance, may alter the spectral signature of an infected plant. To detect 
small reflective differences between healthy and diseased plants, radiometric 
measurements were obtained during the 1982 growing season using a field-portable 
spectroradiometer. The Spectron SE 590 (Spectron Engineering Co., Denver, CO) 
measures reflectance simultaneously at 256 wavelengths (M00-1 lOOnm) , storing it 
on digital tape. Readings were recorded from 2.5 m above the plant canopy, 
using a 10 degree field of view. Black/White IR photographs (35 mm, Wratten 12 
filter) were taken either before or after the spectral scan. Spectral reflec- 
tance ratios were obtained by normalizing readings to a BaSO^ standard reflec- 
tance panel. 

RESULTS AND DISCUSSION 

Aerial photographs obtained during the 1981 growing season exhibited subtle 
changes in reflectivities from bush to bush. Although these changes were 
apparent throughout the growing season, they may represent differences in growth 
stages within and among bushes. Color differences were more obvious on color 
film than CIR. Analysis of the color and black/white IR photographs taken dur- 
ing the 1982 season is in progress and may prove to be a more sensitive tech- 
nique. 

Color and CIR photographs taken 2 m above diseased and healthy grapevines 
suggest reflectance differences may occur early in the growing season. Densi- 
tometric readings of photographs of blueberry bushes demonstrated larger reflec- 
tive differences early (June) in the season for color film and mid-season for 
CIR film (Figs. 2 and 3). 

Generally, diseased bushes exhibited decrease reflectance than healt;., in 
all wavelengths, especially the infrared region. 



Figure 1. 
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FIGURE 4 

FIELD REFLECTIVITIES OF BBSS'/- INFECTED AND HEALTHY CV. JERSEY HIGKBUSH BLUEBERRY BUSHES 
COQPERSVILLE/ MI - 1932 



392.1 497.0 609.7 726.8 849.5 976.0 1104.2 


Wavelength (nm) 
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Spectrophotometric recordings were taken late in the 1982 growing season of 
5 pairs of healthy and diseased plants in fields containing the three virus 
diseases; BBLMV, BBSSV and PRMV. In three of the five pairs observed, BBSSV 
Infected blueberry bushes showed a lower reflectance in the infrared region 
(726.8-1007.7 nm) than their healthy counterparts (Fig. 4). One pair showed the 
reverse situation while the fifth pair showed no difference (Fig. 5). All gra- 
pevines infected with PRMV tested, exhibited decreased reflectance in the 
infrared wavelengths as compared to healthy. The reverse was observed for BBLMV 
Infected blueberry bushes, where four of the five pairs showed increased reflec- 
tance over healthy. One pair showed no difference in reflectance. 

The development of a reliable method for differentiating between diseased 
and healthy plants requires further field and laboratory analysis. Spectropho- 
tometrlc analysis using the field portable unit obtained late in 1982, will be 
made during the entire season in 1983. Both aerial and low level detection will 
continue to be made using various techniques. 

Objectives for 1983 

1. Visual and densitometrlc analysis of aerial and low level photography. 

2. Use of the field portable spectrophotometer at various phenologlcal stages, 

to determine maximum differences between healthy and diseased plants at 
wavelengths ranging from 400 to 1100 nm. 

3. Continued prediction of rate and pattern of spread of these diseases under 


field conditions. 



Spectral Energy (Amplitude) <M 
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FIGURE 5 

FIELD REFLECTIVITIES OF BBSSV- INFECTED AND HEALTHY CV. JERSEY HIGHBUSH BLUEBERRY BUSHES 
COOPERSVILLE, MI - 1982 




CONTRACTUAL ACTIVITIES 


Important Farmlands Inventory 

The Center for Remote Sensing Is oontinuing work under contract with the 
Soli Conservation Service in preparation of Important Farmlands maps for coun- 
ties in Michigan. The mapping involves the delineation of prime soil areas frost 
soil survey information and the identification of unique farmland , water and 
urban built-up areas form aerial photography (U.S. Department of Agriculture, 
Secretary's memorandum number 1827, revised, October 30, 1978). Unique farm- 
lands are lands other than those designated prime that are used for the produc- 
tion of specific high-value food and fiber crops (e.g., tree and bush fruits, 
vineyards and vegetables) . 

The information is compiled onto a 1:50,000 base map of each county and 
area statistics per oategory are determined. The Important Farmland maps are 
being produced under the Land Inventory and Monitoring (LIM) program of the U.S. 
Department of Agriculture. Contracts for the preparation of Important Farmland 
maps of 11 additional counties were undertaken this reporting period, which 
brings the total under contract to 36 counties. 

Changes in Aquatic Vegetation in Saginaw Bay 

During the past decade, concern had been raised over the loss of wetland 
habitats along the shores of the Great Lakes, especlslly the cattail marshes 
which provide food and a cover habitat for waterfowl. Because coastal marshes 
are dynamic environments, changes in the amount and composition of aquatic vege- 
tation may occur frequently in response to changes in water level, site moisture 
conditions, wave action from storms, water chemistry, sedimentation, hunan 
activities, and other processes. Fluctuations in the water level of the Great 


Lakes are a major faetor influencing marsh oonditiona, particularly vegetation 
growth and leas. 

Center for Remote Sensing investigators studied coastal wetland changes and 
docuoented over time the types, distributions, and smounts of aquatic vegetation 
existing within three sites along the shoreline of Saginaw Bay in Lake Huron. 
This work was performed for the East Central Michigan Planning and Development 
Region. 

The project mapped aquatic vegetation within three areas for selected years 
through interpretation of historical aerial photography. Vegetative area by 
category was determined per year and its change in subsequent years. The areas 
were chosen based on the availability of aerial photographic coverage, wetland 
location within Saginaw Bay, composition of aquatic vegetation, and other physi- 
cal factors related to the nature of the marsh. The three areas selected are 
state wildlife areas of Quanicassee, Nayanquing Point and Wildfowl Bay. 

Aquatic vegetation maps were prepared for Quanloassee and Nayanquing Point 
from aerial photography taken in 1950, 1969, 1975, and 1978; and for Wildfowl 
Bay from 19*»9, 1975, and 1978 photos. The selected years of photo coverage pro- 
vide a record of marsh conditions during periods of medlun to high lake levels. 

A vegetation classification system and interpretation criteria were esta- 
blished. Through initial photo analysis, it was determined that four categories 
of aquatic vegetation could be identified on the aerial photographs: Submer- 
gents, cattails, mixed emergents, and sedges and grasses. Cattail areas were 
further classified into three stand density or cover classes (255, 25-755, and 
755), since cover is an important factor in the use of wetland habitats by 
wildlife. 

Eleven aquatic vegetation maps were prepared and then digitized using an 
electronic digitizer. The resulting polygon files were converted to raster 
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(l.e., grid cell) arrays with a 0.1 sort call siza. Thaaa fllaa wara dlaplayad 
on a color monitor for visual analysis and araa statist los wara calculated for 
eaoh map category and listed on a printer. 

The findings of the study Indicated that between 1949 and 1978* Wildfowl 
Bay lost about 401 of Its aquatic vegetation (2364 acres), about twloe as much 
as either Quanlcassee (219) or Nayanqulng Point (179) lost. More Importantly, 
Wildfowl Bay lost a greater percentage of Its hlgh/medlun-density cattail stands 
(359) compared to Quanlcassee (319) or Nayanqulng Point (209). 

These greater losses may have been due In part to differences In wave- 
action protection, since both Quanlcassee and Nayanqulng Point are less exposed 
than Wildfowl Bay. The location of Wildfowl Bay on the eastern side of Saginaw 
Bay subjects it to enhanced wave-action and on-store lake Ice movement under the 
Influence of the prevailing westerly winds. Other factors such as sedimentation 
or water quality differences may also have been Involved. 

A report entitled, "Changes In Aquatic Vegetation In Quanlcassee, Nayanqu- 
lng Point, and Wildfowl Bay," was prepared by Bill Enslln and Ifcayne McIntosh 
for the East Central Michigan Planning and Development Region. 

Digitized Soil Association Map of Michigan 

A soils database for the entire state of Michigan was eompleated this 
reporting period. The database was prepared for Dr. Maurice L. VI tosh of the 
Department of Crop and Soil Science, Michigan State diversity. The base map 
used was the Soil Association Map of Michigan, which was compiled by the 
Cooperative Extension Service and the Agricultural Experiment Station of Michi- 
gan State University and the Soil Conservation Service, United States Department 
of Agriculture. The map base was approximately at a 1 to 1,004,000 scale and 
contained 78 different soil associations. The database resulting from this map 
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contains 633 columns by 733 rows of grid cells. Tht size of ons grid osll is 
ons kllomstsr. A steond vsrisbls, oounty boundaries, was also dlgltlzsd in 
addition to tha soil varlabla. 

This databasa was than usad by Dr. VI tosh and tha Cantar to avaluata tha 
spatial distribution of tha affaotlvanass of soil nitrification inhibitors. Tha 
first stap was to group tha soil typas into flva suitability ratings with 1 
balng a poor rasponsa and 5 balng a good rasponsa to nitrification inhibitors. 
This was dona for tha following tlma parioda: aarly fall, lata fall, aarly 
spring, lata spring and alddrass. Araa totals and peroantagaa wars than gen- 
aratad for thasa groupings, both on a oounty and stata laval basis. Arsas of 
high potential wars than windowed out of tha statewide file and photographed 
from tha color monitor. Tha resulting color enlargements gave good visual docu- 
mentation to supplement tha generated statistics. Tha araa totals and color 
enlargements ware than used by a local ohamleal company to avaluata tha affao- 
tlvanass and market potential of a nitrification inhibitor they ware in tha pro- 
cess of developing and marking. 

Aerial Photography Support for Hybrid-Poplar Research 

Tha Cantar for Remote Sensing acquired aerial photography of hybrid-poplar 
research stands for Packaging Corporation of America, Manistee, Michigan. Tha 
photographs ware evaluated ss part of a larger research effort of tha company to 
access woody biomasa production on a commercial scale. Tha objectives of this 
research are to identify short-rotation intensive hardwood silviculture produc- 
tivity rates from already established large-scale hybrid poplar plantings; iden- 
tify major plantation pasts, diseases and their impacts and control measures 
needed to Insure high biomass productivity in a large-scale operation; evaluate 
the need for weed control in the second growing season; determine the effeet of 
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several herbicides used for post-planting weed control and their effect on tree 
growth in the field; test planting stock and site modification techniques which 
might increase initial establishment success leading to Increased growth and 
survival and present a consolidated economic evaluation of large-scale biomass 
plantations. Light aircraft, 70mm aerial photographs (both color and color- 
infrared) were acquired at four time intervals during the sunmer of 1982. These 
photographs were provided to the various investigators (from Packaging Corpora- 
tion of America, Michigan State University, U.S. Forest Service, Iowa State 
University and Michigan Technological University) for use and evaluation in 
their respective research endeavors. Additionally, Center staff have provided 
limited technical support in the incorporation of the aerial data into the ongo- 
ing research efforts. 

CRIES Project Assistance 

The Center for Remote Sensing assisted in creating a database of Jamaica. 
The database was completed last September and was created in conjunction with 
the Comprehensive Resource Inventory and Evaluation System (CRIES) project. 
Department of Resource Development, Michigan State University. The database 
contained three variables: 13 political divisions (parishes), 23 land use types 
and 240 soil types. The map base was at a 1 to 50,000 scale. The final data- 
base resulting from these maps contained 1220 columns by 480 rows of grid cells. 
The size of one grid cell was 4 hectares (200 x 200 meters). The database was 
initially created at a much finer resolution ( .25 hectares per grid cell) in 
order to capture all of the important information from the base maps. This 
higher resolution data was then aggregated into the final database resolution of 
4 hectares per grid cell so that a more general analysis could be done. The 
high resolution file was retained for future analysis. 
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This aggregated database was analyzed using various techniques. Overlays 
were run to find information on coffee and banana production potential. Many 
different types of area totals were also generated. After the analysis had been 
completed, the database was transfered from our microcomputer to the 
university's mainframe computer. There it was reformated and written out to 
tape in a format compatible with the geographic information system on the main- 
frame computer resident in Jamaica. 

Projects similar to the Jamaica project are in the process of being com- 
pleted for the countries of Kenya and Honduras. In addition to the analysis 
mentioned above, a powerful crosstabulation program has been written for CRIES 
by the Center to further analyze the databases of these new projects. 


EDUCATION AND TRAINING 


Interpretation of Color Infrared Alrphotoa for 
Forest Resource Inventories 


October 28-30, 1980 
June 16-18, 1981 
September 22-24, 1981 
September 27-30, 1982 


Escabana, Michigan 
Escabana, Michigan 
Sault Ste. Marie, Michigan 
Grayling, Michigan 


This three-day shortcourse introduced its participants to the application 
of mediun-scale (1:24,000) color infrared aerial photography to forest cover 
type mapping. In addition to lectures, ample time was allotted for discussion, 
demonstrations and practical exercises. A half-day session was devoted to field 
verification of the interpretation exercise. 

Staff Members William D. Hudson and David P. Lusch were the Instructors. 


Aerial Photography for Natural Resource Management 

March 29 - April 7, 1982 
Kingston, Jamaica 

Dr. David P. Lusch, Center for Remote Sensing, was one of three instructors 
for an 8-day workshop held at the University of the West Indies. Sponsored by 
U.S.A.I.D. in cooperation with Michigan State University, Comprehensive Resource 
Inventory and Evaluation System (CRIES) Project, this training activity 
acquainted the participants with the practical use of aerial photography for 
natural resource management. Fourty hours of classroom instruction (lectures 
and exercises) were supplemented with one day of field work dealing with 
ground-truth collection and interpretation verification. A 217 page study guide 
was prepared which included a set of specially prepared stereograms illustrating 
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the major land cover/use categories encountered in Jamaica. 

Remote Sensing for Land Use Analysis 

Geography 411 or Urban Planning 800, 4 credit hours 
August 2-6, 1982 

Three staff members from the Center for Remote Sensing, W.R. Enslln, W.D. 
Hudson and D.P. Lusch, conducted thl3 5-day, 4-credit course which emphasized 
the use of remotely sensed data to Inventory land cover/use. lhirty-six hours 
of classroom instruction, including none: ous exercises, were supplemented with 4 
hours of field work collecting ground-truth information. As part of the course 
materials prepared for this ola s, the instructors wrote the Photo Interpreta- 
tion Key To Michigan Land Cov er /Use . Thirty-six students attended this course. 

Land Use Mapping From Aerial Photography 
August 10-12, 1982 

This 3-day shortcourse was sponsored by the Center for Remote Sensing, in 
cooperation with the Departments of Geography, Urban Planning and Forestry, the 
Agricultural Experiment Station and the Cooperative Extension Service at Michi- 
gan State University, as well as the Michigan Resource Inventory Program, 
Department of Natural Resources. Emphasis was placed on the practical applica- 
tion of airphoto interpretation for land cover/use inventories with special 
reference to 1:24,000 CIR imagery and the Michigan Current Use Inventory. The 
twenty participants included representatives from the following agencies or 
firms: 


Planning Commissions 10 
Consulting Agencies 3 
Utilities 2 
State Government 2 
Soil Conservation Service 1 
Miscellaneous 2 
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Staff members William R. Enslin, William D. Hudson and David P. Lusch were 
the instructors. 

Michigan State Remote Sensing Newsletter 

The Center produces a quarterly newsletter (the first edition was March 
1981) which is distributed principally in-state. The purpose of the newsletter 
is to keep state, regional and local decision-makers informed about new remote 
sensing capabilities and applications and to reach new potential users. 

SPECTRUM 


SPECTRUM is intended to keep Center investigators and others interested in 
remote sensing and its applications up-to-date on current news in the field. In 
SPECTRUM we attempt to brief readers on such things as recent developments in 
the field (principally excerpts from newsletters), upcoming meetings, workshops 
and seminars; staff travel; planned flights, etc. SPECTRUM is not released on a 
planned schedule— only as it is timely. 



PUBLICATIONS 


of the Michigan State University, Center for Remote Sensing 


1. Applicability of Satellite Freeze Forecasting and Cold Climate Mapping to 

the Other Parts of the United States . Center for Remote Sensing and 
Department of Entomology, Michigan State University, November 1981. 
54p. 

2. "Automatic Digital Image Registration." A. Goshtasby, A.K. Jain and W. R, 

Enslin. Proceedings 8th International Symposium on Machine Processing 
of Remotely Sensed Data, pp. 347-352. Purdue University, West Lafay- 
ette, Indiana, June 1982. 

3. Interpretation of Color Infrared Airphotos for Forest Resource Inventories . 

William D. Hudson and David P. Lusch. Center for Remote Sensing, 

Michigan State University. July 1982. 55p. 

4. Photo Interpretation Key to Michigan Land Cover / Use . William R. Enslin, 

William D. Hudson and David P. Lusch. Center for Remote Sensing, 

Michigan State University. August 1982. 71p. 

5. "Assessment of Modified Surface Temperatures and Solar Reflectance Using 

Meteorological Satellite and Aircraft Data." J. Bartholic, S. Gage, 
A. Goshtasby, C. Mason. Presented at the Symposium on Study of Land 
Transformation Processes from Space and Ground Observations, Ottawa, 
Canada, August 1982. 

6. Changes in Aquatic Vegetation in Quanlcassee , Nayanqulng Point and Wildfowl 

Bay . Prepared for the East Central Michigan Planning and Development 
Region. Bill Enslin and Dwayne McIntosh. Center for Remote Sensing, 
Michigan State University. October 1982. 

7. "An Evaluation of Digital Landsat Classification Procedures for Land Use 

Inventory in Michigan." Richard Hill-Rowley. Ph.D. dissertation, 
Michigan State University, 1982. 258p. 

8. Visual Interpretation of Landsat Imagery for the Identification of Coni- 

ferous Forest Types in Michigan. Kathryn L. Franklin, William D. 
Hudson, Carl W. Ramm, Center for Remote Sensing and Department of 
Forestry, Michigan State University. Submitted for publication to 
Michigan State University Agricultural Experiment Station. January 
1983. 

9. A Methodology for Constructing Vegetation/Physiographic Profiles Using 

Remote Sensing and Collateral Materials. William D. Hudson. Center 
for Remote Sensing, Michigan State University. In preparation. 
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10. Identifying forest Types end Tree Species from Color Infrared Aerial Photo- 
graphs. William D. Hudson. Center for Remote Sensing, Michigan State 
University. In preparation. 


ASSESSMENT OF MODIFIED SURFACE TEMFERATURES AND SOLAR REFLECTANCE USING . 
METEOROLOGICAL SATELLITE AND AIRCRAFT DATA 

J. Bartholle*, S. Cage**, A. Coshtaaby***, C. Mason*“* 

■Agricultural Experiment Station, ■•Oapartoant of Entomology, ■■•Cantar for Re- 
aota Scnalng 

Michigan Scatc Unlvaralty, Eaat Lanalng, Michigan 4882* 

■••■Department of Atmospheric and Ocaanlc Sclanca 
Unlvaralty of Michigan, Ann Arbor, Michigan 48104 


ABSTRACT 

Change* In aurfaca caaparatura raaultlng fron tha actlvltlaa of hii era avaluatad ualng 
oataorologlcal aatalllta (NOAA and HCMM) and alrcrafc data. Study altca vara locatad In 
Florida and Hlchlgan. Tharoal data abound that day aurfaca caaparatura* over large araaa 
could be lncraaaad by 10-iS*C by aodlf Icatloos raaultlng froa agricultural practice*. 

Change* In reflected aolar radiation a* a function of agricultural practice* uaro detectable 
using HOOf data. 1/ 

INTRODUCTION 

Through tins, nan la gradually Modifying oorc of the earth's surface to optlnlz* conditions 

for habitation. The changes era largely for crop production, anlnal grazing and utilization 

of blonat* for cooking and heat. In thi* procesa forest* ora being converted to pastures 
and farmland*, grassland* *ro being grazed by dooastlc animals, nau specie* of plants are 
being Introduced and major portion* of the surface are being drained or Irrigated. 

Thl* process has evolved slowly over centuries, but has speeded up dramatically during thi* 
century. Increasing oodlflcatlon of the surface for food production has been an essential 

part of increasing the earth's carrying capacity from two billion at the turn of the century 

to over four billion at the present doc. Further, thete changes will be occurring at an In- 
creasingly rapid rate aa the projected population of the earth doubles by the year 2020 (!]. 

Relatively few studies have assessed the magnitude .of these changes or their Impact on sur- 
face boundary conditions (2,4). Many characrerlstlcs of the surface, Including temperature, 
radiation and surface roughness, could be altered significantly. This paper shows that in- 
formation available fron the National Oceanic and Atmospheric Administration (NOAA) and NASA's 
Neat Capacity Mapping Mission (HOW) satellites and airplane scanner data can help assess the 
Impacts on surface temperature and radiation of man's modification of our planet. 

METMODOLOCY 

Several recently launched meteorological satellites have collected valuable data on re- 
flected radiation In the .5 to 1.1 micron range. This band has the Information needed for 
good estimates of aolar reflectance from different surfaces (5). Also, data froa scanners 
sensing In the 10. 5 to 12. 5 micron range are available on several satellites which can be 
used to determine surface temperature. Accuracy and resolution of these systems Is suffi- 
cient to characterize vegetative types and to provide answers about the magnitude of changes 
that might be ezpected with deforestation or other significant change* In existing vegetation. 

NOAA eatelllte data provide* on* kilometer resolution In the thermal (10.) to 12.5 microns). 
Th« satellite orblcs over at approximately 1000 and 2200. The sun-synchronous HCKH satellite 
provides high resolution thermal data with a .6 km by .6 km reaolutlon at nadir. In tha 10.5 
to 12.5 micron range. Also, the concurrent visible channel 0.5 to 1.1 microns with a dynamic 
rang* of approximately 100 percent albedo and resolution of .5 km by .5 km at nadir poten- 
tially will provide the essential data required to meet the objective of characterizing aolar 
reflectanca (6). 

For tha NOAA data, tranaparenelea were used, digitized end analyzed on the Image 100 at the 
Kennedy Space Center. The HCMM data was obtained In both pictorial and digital forms. Ana- 
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lysis was toa at the Cantar far lamia Sensing inf tha lamia Processing Laboratory at HUM- 
gaa Stair University. 

Tha digital NOW raflaetad radiation data warm uaad to find tha ralatlon between vaiatatlon 
types and raflaetad valuaa. Tha raflactanea laaga waa thraaholdad at dlffatant values aad 
tha regions shoving tha aaaa raflactlvlty vara laolatad for aach claaalf Icatlon. Tha iae- 
latad regions war a coopered with an alraady claaalf lad Landaat image. 

A Oaadalua 01-1250 acannar with 8-1* microns thanaal channal waa flown to obtala hl|h raso- 
lutlon inputs required to aaolat In tha analysis. Thla data waa acquired on analog tape 
than dlgltitad and analytad on tha Image 100 at tha Kannady Spaca Cantar. 

Tha atudy araaa vara mainly In Florida and Michigan. Tha NOAA aatalllta and alreraft data 
wara ohtaload for Florida and HCMM data waa uaad for analyale of tha Michigan taae alta. 

RESULTS 

NOAA thermal aatalllta lwgary for April 28. 1978 at about 1000 (Fig. 1) ahowa darker araaa 
aa warmer and lighter araaa aa cooler. The agricultural araaa In tha Klaalaalppl Valley 
and tha a out here portlooa of tha United Staten Including Ceorgla, Alabama and South Carolina, 
plua major araaa of Florida, are ohown to ba markedly warmer (darker) than the adjaceot 
naturally vegetated araaa. Separate analyale lndlcatea that many of the darker araaa are 
about 10-12*C warmer than tha adjacent natural araaa. 

Mora detailed emaalnatlon of tha Florida penlnaula la poaalble In Fig. 2 (note that warmer 
araaa appear lighter, tha raweraa of Fig. 1). In April, many flelde In north central Florida 
are tilled and cropa do not yet cover a elgnlfleant portion of tha aurfaca. Tha a a agricul- 
tural araaa are In long, broad atrlpa generally lying In a north-south direction and arc much 
warmer. Soma of tha greateat thermal contraata In eouthern Florida occur between tha water 
cooeervatlon raaervaa In tha Evergladeo and tha adjacent agricultural araaa that have bean 
drained. Tha adjacent agricultural araaa are to the north and juat eouth of lake Okeechobee. 
Seth drainage and agricultural, practlcaa are affecting tha aurfaca condltlona, caualng 
taaparaturea to ba elgnlflcantly warmer than natural araaa. 

Image dleplaya af aircraft thermal acannar data are ehovn in Flgvrea JA and 31. Two dls- 
almllar araaa wara choeen for atudy in tha Taylor Creak watcrahed located north of Lake 
Okeechobee, Florida. Tha firat area, ehovn in Fig. 3A, waa primarily Improved graaa 
paetures uaad for gracing dairy cattle. Tha aacond area (Fig. 18) lncludea paatura, march 
and a large cltrwa grove. The thermal data of Fig. 3A wara obtained at 1222-1223 1ST on 
April 28, 1978. Wlndapeed waa 700 cm/aec, air temperature waa 28.1*C and tha relative 
humidity waa 39 percent (dew point 11.1*C). Thermal data of Fig. IS were obtained at 1*32- 
1*33 EST on April 28, 1978 and meteorological condltlona ware dollar to thoaa measured an 
April 28. Iheae aurfaca meteorological meaaureaenta wara made at a height of 10m. A uniform 
general rain of about 3cm had fallen on tha area two daye prior to the fllghta. Soil oola- 
tura condltlana were good for growth (3). 

U.S. Highway 4*1 runs N-S through tha Image In Fig. 3A. To tha aaat of tha highway tha 
affecta of high dairy cattle population denalty are easily aeen. In the irregularly- 
a ha pad araaa heavy cattle traffic had eliminated planta. A alsture of aand and partially 
decompoeed manure made up tha aurfaca cover. The aurfaca temperatures In this area wara 
greater than *2*C. Similarly, Field 1, which had apparently bean cut for hay and than 
Srated heavily prior to April 28, ahowed little aprlng growth and temperaturea wara in 
emcaea of *2*C. Field 2 of tha aame apeclea, primarily Fangola ( Dltltavla decumbaue ) had 
bean cut For hay but had not bean graced. It waa in a healthy growing condition and ranged 
from 38-*2*C. Field 3 waa rather heterogeneoua, containing patchae of graaa interapereed 
with broadlaavad weeds and thinner atanda of graaa and ranged in temperature from 2S-38*C. 
Evapatranaplratloa (ET) and heat flue valuaa for Field 1 were 20 and 30 w/ea‘, raapectivaly, 
and far Field 3 with higher ET were *0 and 10 aw/cm 2 (3). Vet araaa with woody ahruba aa 
tha primary pound cover had the cooleat temperaturea In the scene with a range of 22-2i*C. 

la Fig. 38 the marsh vegetation was coolest (about 21-2**C) and the citrus ranged from 
2*-28*C. Feature 1 was poorly managed and had temperaturea from 31-38'C while Feature 2 
waa not ovar-gratad and ranged in teaparature from 2*-31"C. 

Fig. *A shows the MOM thermal Infrared data- over the northern lower penlnaula of Michigan. 
This imagery is from September 28, 1979. Olgltal values showed that the agricultural areas 
ware weraast when compared to ether areas. Surface redlant temperature for agricultural 
areas average 20*C while the average radiant temperature for water waa 8*C. The surface 
temperatures for other surface types fall between theee two valuea. Average tenperaturas of 
surface types acquired from thla scene are shown In Table 1. The distinct differences la 
temperatucea between the aurfaca types ahowa tha relationship of aurfaca temperature and 
surface types. 
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TAIL! 1 Average Ttiptnturi for Dlfforont Surfoe* Tjrpu oa OkMlMf froa MOM 
Data oa September 24, 1129 far tha Northara fart lan af the Lever 
Panlnaula af Michigan 


Surface Type 


Water 

••c 

Swamp 

13 *C 

faraot 

17*C 

Agricultural Araaa 

20*C 


Fig. 41 ia tha DOM raflactanca laago af tha aaaa acana. Again, digital valuaa ehoved tha 
dependency af raflactanca ta tha aurfaca typai. Experiment a ehoved that It la a van pa a a lb la 
ta verify aaaa vegetation typee aolely by thalr raflactanca valuaa. In fig. II raflactanca 
valuaa batvaan .0? and .09 ara laalatad In black, which ahav tha cealfereue arena vary vail 
(aaa tha araaa painted to by the errowa). This vaa verified by referencing tha Laadaat 
claaalfled lnage of Auguac 1, 191} (1], fig. AC. Thle ahawa that the MOM aatalllta raflac- 
tanca data can be uaed ta characterlaa thla earth aurfaca paraaatar. 

DISCUSSION 

Fran tha analyala af tha aatalllta data. It la dear that aan'a actlvltlee have algnlflcantly 
alterad aaaelve portlona af tha aarth'a aurfaca. Theac changaa have, In turn, modified tha 
planat'o boundary temperatvree. Tha change la not- only largo In area, but alaa In magnitude, 
fraauantly l<b»12*C varnar tenpereturee vara found In agricultural arena when coopered ta tha. 
natural aurfaca (fig. 3 and A). Tha a a nadlflcatlana could algnlflcantly lnpact tha tharnal 
radiation leaving tha aarth'a aurfaca and the repartltlenlng af energy Inca aanalbla and 
latent flunaa. 

further, tha alrerafc tharnal data vividly ahav (fig. M and )•) that even vlthln the agri- 
cultural araaa, aa vagctatlen types ara changed or overgrating eccure, there can ba further 
algnlflcant changaa In tha tharnal reglaea. Tha a a acence ah owed ovargratad araaa vara 
greater r1)aa.42*C while wall-watered nondornant paaturaa vara appreeinataly 30*C. far thaaa 
tanparatura dlffarancaa to occur, najar dlfferencea In energy going Into avapecranaplretlen 
bui t anlat. Thue, tha flunaa af vapor and heat ara canaldarably different depending an haw 
the 'aurfaca la managed. Thaaa aurfaca dlffarancaa will ultimately Impact tha hydreleglcal 
balance. 

Tha lapacta af thaaa changaa ovar tha caapaelta area of tha Florida panlnaula could poten- 
tially have algnlflcant lmpacta an local vaathar. Soma evidence of thin modification la 
vlalbla In fig. 3. lake Okeechobee, which la a man-made lake, la ahevo ta ba modifying the 
cloud pattern ovar tha eovthern portion af the Florida panlnaula In thla figure. 

Nadlflcatlana In a mere temperate area arc dearly ehovn by tha example uaing MOM data far 
tha panlnaula of Michigan. Thla analyala alaa chewed that algnlflcant changaa in raflactanca 
could occur with changaa In vegetation. Madlflcatlona In thla reflectant ceaponent could 
have long term ecological lapllcatlena alnca aalar radiation la tha main driving farce far 
taaparatura Inereaac, hydrological changaa and biological pracaaaaa. 

CONCLUSION 

Nadlflcatlana at aurfaca taaparatura and radiation vara clearly ebaervabla from analyala af 
aatalllta data. Both tha epatlal dictrlbutlen and magnitude of thaaa changaa can potentially 
ba monitored uaing aatalllta data. Thue, coma preliminary cenelderetlona for the habitability 
af araaa of tha north can atart to ba developed. Of greataat concern, however, la that tha 
papulation la expected to lncraaaa by an additional SOS by the turn of the century, aa tha 
proc'iaa of change exanlned In thla paper can ba expected to accelerate In Che yearn ahead. 
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Fig. 3 Iftagaa of aircraft thermal acanner data. 
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Fig. ! Ther&al *ccne for the eaatern United Staten from NOAA satellite for 28 April '78 at 
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ABSTfUCT 


This paper introduces a general 
procedure for automatic registration of two 
images which may have translational, 

rotational, and scaling differences. This 
roccdure involves 1) segmentation of the 
mages. 2) isolation of dominant objects 
from the images, 3) determination of 

corresponding objects in the two images, 
and 4) estimation of transformation 
parameters using the center of gravities of 
objects as control points. An example is 
given which uses this technique to register 
two images which have translational, 

rctationjl, and scaling differences. 


I. INTRODUCTION 


Given two images of the same scene .n 
coordinate spaces (x.y) and fx'.y'), iiru.ge 
registration is the determination of 
transformation functions fl and f2 such 
that given the coordinates of a point in 
one of the images, we can compute the 
coordinates of the same point in the other 
image by 

x • f 1 ( x * ,y* ) 
y • f 2 (x * ,y’ ). 

The parameters of transformation 

functions fl and (2 are estimated using a 
set of corresponding control points from 
the two images. Selection of control 
points by hand is often time consuming and 
is susceptible to systematic errors. It is 
desirable to make the control point 

selection process automatic so that the 
whole registration process can be carried 
cut automatically. 


An autosatic technique for selection of 
control .points in the first image is 
available * . This involves selection of 

windows which 1) contain a large number of 
high gradient edges, 2) contain a large 
number of connected edges, 3) are unique 
with respect to neighboring windows, and 
4) are well dispersed in the image. The 
upper left hand corner of the window 
satisfying these properties is taken as a 
control point. In the second image, **•. 
control points are obtained by carrying 
a search for the position of best r * l . f .< 
each window from the first image J,i >*:n' 
either thg sequential similarity ;oh 

algorithm? or the cross-con ..ution 
technique 1 . 

Selection of control points 

automatically in this manner has shown to 
be satisfactory as long as the two images 
have only translational differences. for 
imuges with rotational. and scaling 
differences, the search process becomes 
inaccurate and unreliable. In sueh 
situations, image points which are 
routinely used as control points include 
intersections of lines (intersections of 
roads in aerial or satellite images) and 
positions where lines join (like positions 
where rivers join). 

In this paper we introduce a technique 
for automatically finding corresponding 
control points in two images which may have 
translational, rotational, and scaling 
differences. The technique uses center of 
gravities of objects as control points. 
Because the coordinates of the center of 
gravity of a shape is the average of the 
coordinates of the pixels on its boundary, 
if any error has been made in extracting 
the boundary, that error is averaged over 
the whole boundary and so the effect on the 
center of gravity is small. Another 
beneficial property of the center of 

ravity as control point is that its 
ocation can be determined upto subpixel 
accuracy while the traditional control 
points (like intersection of roads) can 
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take only discrete values. Our 
registration algorithm consists cf 
following stops. 1) Segmentation of the 
images. 2) isolation of dominant objects 
from the images, 3) determination of 
corresponding objects in the two images, 
and 4) estimation of transformation 
parameters using the center of gravities of 
objects as control points. 


II. IKA6S SEGMENTAT 1 ON 


There is no single image segmentation 
algorithm which can reliably segment an 
arbitrary image. Depending on the type of 
imagery, different techniques have been 
developed by different authors. Some of 
the techniques like the recursive region 
splitting technique’., the semantic region 
growing technique , 7 or the vgion 
classification technique work especially 
good on multispeetral aerial or satellite 
images. 

For HCMM satellite images which are of 
interest to us, a simple segmentation 
technique based on gray level thresholding 
appears to give satisfactory results. A 
number of techniques have been proposed to 
obtain an appropriate value of the 
threshold. The threshold value can be 
selected by minimizing the global average 
contrast of the image 0 . In another 
method , the threshold value is determined 
by minimizing the sum of the squares of the 
differences of the original image and. the 
thresholded image pizel by pixel. Katz has 
determined the threshold value using the 
edge information present in the image. We 
have chosen the Katz algorithm for the 
segmentation of HCKM images because of its 
speed and performance on this type of 
image. The algorithm consists o' the 
following four steps. 

1. Compute the gradient of the original 
image, call it image C. 

2. rind the high gradient pixels by 
thresholding image C. Replace all 
values above the threshold by 1 and all 
other values by 0. Call the new image 
N. 

3. Multiply H by the original image, call 
the new image M. Note that M will 
contain only high gradient pixels of the 
original image. 

4. Compute the average gray value of image 
H (nonzero pixels only). This will be 
the required threshold value. 

Note that the algorithm requires a 
user-specified parameter in step 2. 
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III. FINDING CORRESPONDING REGIONS 


Once the two images are segmented, ve 
isolate regions with closed boundaries in 
each one of the images. A region which 
touches the boundary of the image is not 
isolated, rurther, very small regions are 
discarded. The task is now to determine 
the correspondence between the two sets of 
regions. 

Regions in the two images could be 
matched based on shape alone. A number of 
shape matching procedures are available in 
the literature, fourier descriptors, and 
invariant moments * 1 have been extensively 
used for shape matching. The distribution 
of chord lengths joining the center of 
gravities of regions to its boundary polsts 
has been proposed as a measure of shape . 

Price and Reddy 14, ^extracted a number 
of features from each region including 
roundness (perimeter /4ir area ) . length to 
width ratio, color, intensity, location, 
relative position to mateh the regions is 
the two images. Their matching technique 
seems to best fit our problem and we have 
used it with some modifications in the 
example given in section V. 


IV. TRANSFORMATION FUNCTIONS 

Selection of the right transformation 
function is another important factor in the 
registration of digital images. The best 
transformation function for registering two 
images which have only translational 
differences is very simple and has only two 
unknown paraeteters. Applying a 

transformation function with more 
parameters makes the registration process 
costly and probably less accurate. On the 
other hand, images of the same scene taken 
at different angles require a 
transformation function with I unknown 
parameters (the projective transformation) 
»i.d a transformation function with less 
parameters cannot register them accurately. 

Usually it is known (or assumed) that 
the two images have only specific 
differences and can be registered using an 
appropriate transformation function. There 
are cases theugh, where no a-priori 
information is available about the images 
and we have to select the best 
tra: sformation function to register thea. 
Using the control points, it is possible te 
find some knowledge about the images and 
then select the appropriate transformation 
function. For example, 

1. Using two pairs of correspond* ng control 
points, it is possible t find two 
corresponding line segmentss in the tve 
images. If the ratio of corresponding 
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line segments is not equal to 1, it 
shoes that the tuo images have scaling 
differences. 

2. If the ratio of corresponding line 
segments are not the same for different 
line pairs, it shoes that one of the 
images is g metrically distorted with 
respect to the other. 

3. If the angle between corresponding line 
pairs in the two images are different, 
it shows that the images have not been 
obtained from the same angle. 

The above information helps in the 
selection of the right transformation 
function. 

To be able to estimate parameters of a 
transformation function with n unknowns we 
need at least n/2 corresponding control 
points (which are not eolinear). Usually 
more control points are used and parameters 
are estimated by minimizing the mean-square 
error . 


V. Ml EXAMPLE 


To show how the proposed technique 
works, we have taken a 240x240 subimage of 
the day-visible image obtained by the HCMM 
satellite on 26 September 1?' , 9 from an area 
over Michigan(scene id:A-A0518-18110-l ) . 
The digital image acquired by the satellite 
(figure 1) was used as the first image. 
The second image was the digitized image of 
the print of the same scene provided by the 
National Space Science Data Center. We 
arbitrarily digitized the print so that the 
second image is translated, rotated, and 
!u.s scaling differences with respect to the 
:n st image, as shown in figure 2. The 
digitization was done using the Spatial 
Data System' Vidicon in the Patter-' 
Recognition 6 Image Processing Laboratory 
of the Computer Science Department, 
Michigan State University. 


A. SEGMENTATION 

The Katz technique's used to segment 
the images. The gradient images were 
clipped at 981 of the gradient histogram 
area. Thus average intensity of the 2% of 
the highest-gradient pixels was used as the 
threshold value for segmenting the gray 
level images. We found that other 
day-visible images of HCMM can be segmented 
satisfactorily using the same parameter 
value (98%). Figures 3 and 4 contain the 
segmented images of figures 1 and 2, 
respectively. 


R . CONTROL POINTS 


Ten control points were used to register 
the two images. To obtain the control 
points, we have used a technique similar to 
the one proposed by Price and Reddy* to 
find corresponding .regions in the two 
images. 

Ooject size (perimeter) and shape 
(roundness-perimeter / 4*area) were used to 
find the two most similar object pairs in 
the two images. Let's call them P and Q. 
Now relative distance (distance of the 
center of gravity of an object to the 
center of gravity of P/distance between 
center of gravities of P and Q), and 
relative position (angle beteen the line 
connecting the center of gravity of an 
object to the center of gravity of P, and 
the line connecting the center of gravities 
of P and C> are used to match the rest of 
the objects. Since we have two kinds of 
objects (bright objects in dark background 
and dark objects in bright background) , 
the match is carried out only between 
objects of the same kind. Once the 
rorrespcnuing objects in the two images are 
determined, their boundaries are extracted 
(see figures 5 and 6), and their center of 
gravities are computed which are the 
control points. 


C. TRANSFORMATION FUNCTIONS 

The two images which we are registering 
have been derived from a single image by 
artificially introducing translational, 
rotational, and scale change. Therefore, 
one would suspect that the transformation 
of Cartesian coordinate systems would be 
able to bring the images into registration. 
Using the control points, we computed the 
ratio of corresponding line segments and 
found that the ratios change slightly when 
changing the line pairs. This shows that 
one of the images has small geometric 
distortion with respect to the other. Such 
distortion might be due to the digitizer 
lens or shrinking of the print's paper. 
The transformation function that can 
register images with geometric distortions 
is the polynomial mapping function. We 
used a second order polynomial mappino 
function given by, 

x ' ■ a0*alx*a2y*o3xy*a4x 2 *a5y 2 

y ' ■ b0*Llx-b2y*b3xy»b4x 2 *b5y 2 

Figure ? shows a resampled image of figure 
2 using the above polynomial mapping 
function with the nearest neighbor 
technique. We also used the transf ormat ior. 
of Cartesian coordinate systems and the 
affine transformation. However, these 
transformation functions did not give 
satisfactory results. 
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The accuracy of the whole registration 
process to a large extent depends upon the 
accuracy of control points. In order to 
check the accuracy of the control points we 
segmented the images with different 
threshold values (from 90% to 99.5%) and 
coaiputed the center of gravities of the 
objects for each threshold value. The 
maximum shift in the center of gravity was 

0.7 pixels and the average shift was 0.4 
pixels. Therefore, the center of gravities 
appear to be stable over a reasonable range 
of threshold values used in the 
segmentation process. 


VI. CONCLUSION 


An automatic technique for registration 
of images that may have translational, 
rotational, and scaling differences has 
been presented. This technique can also be 
applied to iatages with small shearing and 
geometric 'distortions (such that in the 
domain of an object the shearing- and 
geometric distortions is negligible). 
Images with distortions due to earth 
rotation, earth curvature, and scanner 
nonlinearity can be registered using this 
technique. 

More research is required to see the 
feasibility of this technique in 
registering images obtained by different 
satellites and/or different sensors. 
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figure 1. of f igure 2. 



figure 5. Dominant objects of Figure 6. Dominant onjects of 

image of figure 3. image of figure 1 4. 
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Figure 7. Resampling of image of 
figure 2 to register with image of 
figure 1. 
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